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ABSTRACT 

The number of <5 rays with energies between 50 and 150 keV that are 
produced by heavy nuclei in emulsions is calculated. The Z 2 dependence pre- 
dicted by the first Bom approximation is corrected by a direct calculation of 
the Mott exact phase^shift scattering cross section. Comparisons are made with 
corrections predicted by the second Bora approximation. When the phase-shift 
results are applied to the problem of charge identification, corrections of up to 
4 units of charge for 1.457-GeV/nucleon nuclei with charge Z = 75 are found. 



INTRODUCTION 


Fowler et al. have recently published results on the detection of heavy 
nuclei in cosmic rays with charge Z > 64 with an estimated accuracy of 2 per- 
cent at the geomagnetic cut-off energy of the experiment [i], The analysis of 
the delta rays produced in the emulsion was based upon the first Born approxi- 
mation. Corrections to the Z 2 dependence predicted by the first Born approxi- 
mation were made on the basis of published results [2 ] of the ratio of the Mott- 
to-Rutherford scattering cross section at the minimum observed 6-ray energy. 

Since the first Bora approximation condition Z/137/3 < < i is not ful- 
filled for large values of Z, Semikoz suggested that a reduction in the reported 
charge estimate by Fowler was necessary, [3]. He based his results on a 
calculation utilizing the second Bora approximation. 

The purpose of this paper is to present an investigation of the determi- 
nation of the charge of heavy nuclei using the Mott phase -shift calculation for 
the scattering cross section. After this cross section is calculated, a numeri- 
cal integration is performed over the delta -ray energy range to obtain the 
number of delta rays produced per unit path length of the heavy nucleus. These 
results are then compared to, those obtained using the first and second Born 
approximation . 


THEORY 

Because only those delta rays with energies much greater than the bind- 
ing energies of the emulsion atomic electrons are used for charge identification, 
one may assume that the atomic electrons are initially free and at rest. Thus, 
by a transformation to the frame in which the electron is initially at rest, the 
cross section for delta fay production da/dT may be obtained from the Mott 
cross section da/d ft normally used for the scattering of an electron from a 
point nucleus. One then finds the cross sections transform as 

da _ 4 7T da 

dT “ T dfi * 

max 

where T, the energy transfer to the electron, and the maximum energy 

transferable to an electron, are related to the electron-scattering angle in the 

cross section differential with respect to solid angle by the relation 

sin 2 0/2 = T/T . If the nucleus is considered infinitely heavy in comparison 
max 

with the electron, then T w 2 m c 2 R 2 (1 - fl 2 )” 1 . 

max e 



In the first Bom approximation, the cross section is 



_ 2tt r n 2 Z 2 / , 
p 2 t 2 y 1 


- p 2 T/T 

max 



( 1 ) 


where r 0 is the classical electron radius and Ze is the charge of the incident 
nucleus. 1 


In the second Bora approximation, the result is 



where ot = 1/137 is the fine -structure constant. 

Therefore, one has a Z 2 charge dependence in the first case with an 
additional term proportional to Z 3 in the latter case. However, since the 
expansion in (a Z) converges slowly, one should consider the higher-order 
terms . 


An alternative approach is to use a phase-shift calculation. The cross 
section thus obtained for the case of a point coulomb potential is [4]; 

(«)„- ^ [ |F ' |JT ' ,+ |G|2 ( T max- T )"] • <3) 

where 

F' = i q F/y , F = F 0 + F t , G = G 0 + Gi 
y = ( 1 - /3 2 ) -1 ^ 2 , q = aZ/p 

F • = 2 ru+jq) ex P [‘ q 111 ( T/T max)] 

1. Unless otherwise indicated, m c 2 units are used throughout this 

6 

paper, where m. is the rest mass of the electron and -ft = m = c = 1. 
e e 
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G » = - 1< l[( T ma X - T ) /T ] F » 

00 

Fi = \ Z A P (p) 

■ f=0 
00 

7 Z B P (p) 

2 f =0 

(-i) i [f 2 ^ - (i +1)^ +1 ] 

ejq? (-iff p £ ) r(p i - i q) 
(pf+iq) F(p i + iq) 


T/T implies that cos 0=1-2 T/T « 
max max 

Two transformations were made on the F t and G t series to obtain more 
rapid convergence . First, to improve the convergence for small energy trans- 
fers, the reduced series employed by Sherman [4] with m = 3 was used. This 
transformation is 

(i -p) m F 1 = ^ Z A (m) P i (p) , 

1=0 

where 

(m) (m-1) 1+ 1 (m-1) i . (m-1) 

i A i "2£+3 A , “ A f-1 

i + i 2i - 1 

and similarly for Gj. 


n = exp, (-i 7ri) r(i - i q) 
l (i + i q) T(i + i q) 

p, = Ji 2 - (a Z) 2 


and 


p = 1 - 2 T/T 


max 


which is found by noting that sin 2 Q/% = 
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Next the Euler transformation, which in our case results in a more 
rapidly converging series, was made to improve the convergence for large 
energy transfers. This transformation is 


Z(-i ) 1 a,- 


i=0 


' £ o 2 k+1 „=0 » ' 


or, in our case, 


ki 


A< „ m) V"> 


. 00 , k 

p _ i y _l_ y . 

1 2 , « „k+i u n (k-n)Ini n v m 

k=0 2 n=0 (i-ji) 


. 00 . k P (n) 

_ _i_ y 1 y kl, n n 

1 2 , „k+i (k-n)I nl . m 

k=0 2 n=0 (i-p) 


The number of delta rays produced with kinetic energies between T t and 
T 2 per centimeter of the cosmic -ray path is found by calculating 



where 17 is the 'number of atomic electrons per cubic centimeter of the emul- 
sion. In the case of the first and second Born approximation, this calculation 
can be done analytically. However, the phase-shift result must be obtained 
numerically. 


Following Semikoz, after the number of electrons produced is found 
from Equation (4) , the comparison with the first Born approximation is made 
by computing [3]: 


AN= (N. - N )/N. , (5) 

1 JD 1 1 

where i implies either the second Born or phase-shift calculation and N is 

B1 

the number of delta rays obtained by using the first Born approximation. 
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An estimate of the charge correction for results based upon the first 
Born approximation can then be found simply as 

Z R =Z (i-AN/2) , (6) 

or |AZ| = | Z - Z | , where Z is the charge found by using the first Born 
It 

approximation. 


CALCULATION PROCEDURE 

To calculate the phase-shift formula, a computer code was developed. 
The approach that was used followed that of Sherman in his calculations of 
electron scattering from point nuclei and his results were duplicated [4 ] . In 
addition, the small-angle approximation of Bartlett and Watson was utilized 
when very small angles were required by the specific case being run [5 ] . 

Double precision arithmetic (a resolution to 18 significant digits for 
the UNIVAC 1108 system) was used throughout the calculation. Table 1 
illustrates a sample test that was used to check convergence. A Newton-Coates 
technique was used to perform the numerical integration with a resulting error 
of less that 1 percent. The cross section values calculated are estimated to be 
accurate to 3 significant figures. Round-off error was estimated by comparing 
the numerically summed series on the left-hand side of equation (7) to the 
analytically equivalent one on the right -hand side [6], 

00 

Z (2i + 1) P (n) x*= (1 -x 2 ) (1 - 2ji + x 2 )” 3/2 , (7) 

i =0 

with x = e~°* 01 and i = 0 to £ = 175 . Agreement was found to be within four 
significant figures . 

Stieltjes' continued fraction was used to compute the gamma functions 
of complex argument [7]. The Legendre functions were computed by using a 
recursion formula with P 0 (ju) =1 and P t (,/u) = ft. 

Since a very large number of terms are required at small angles, 
Bartlett and Watson [5 1 derived an approximation to evaluate the Mott phase- 
shift formula as 0 -*■(). This expression transformed to the rest frame of the 
electron is given by: 
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( 8 ) 



2tt r» 2 Z 2 

pZ -p2 


[ 1 + lr “ z (:r ) 

U • m nv ' 


max 


1/2 


cos 


y 



where 


iy_ r(l/2 - i g) r(l+i q) 

r(l/2 + iq) T(l-iq) 

iy 

or cos 1 y = Re (e ) . 

Whenever T/T corresponded to a small angle, the largest value of 
max 

T for which 



|'d£\ 1 

/(&£.) 

L\ d T/gA 

^PS-V 



was found. For all smaller values of T, the small -angle approximation was 
used, whereas for all larger values the regular phase-shift result was used. 
The two cross sections are shown in Figure 1 for the cases Z = 26 and Z = 104. 


RESULTS 

Figure 2 shows the ratios of the phase-shift cross section to the first 
and second Born approximations for Z = 92 with ji = 0.92. It may be noted from 
the figure that the second Born is higher than the phase-shift result throughout 
this range of energy transfer. This is not the case for an iron nucleus, as 
shown in Figure 3 . 

Figure 4 is the plot of the ratio of the second Born to the Rutherford 
cross section and the phase-shift calculation to the Rutherford cross section as 
a function of Z for nuclei with a kinetic energy of 3.0 and 1.457 BeV/nucleon. 2 
The energy transfer was assumed to be 50 keV. In the figure the ratio for the 
phase-shift calculation reaches a maximum at around Z= 60 for 0 = 0.92 and 
at Z = 75 for (3 = 0.9712, and then decreases at the higher Z nuclei. However, 
the ratio for the second Born approximation increases in proportion to Z as 
expected. Therefore, in the energy ranges that are considered here, the 


2. To convert electron volts to SI units in joules, multiply by 
1.60210 x 10“ 19 . 
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higher -order terms in (aZ) become important whenever (aZ) is greater than 

about (50/137). The third Bom approximation, i. e., retention of terms 

through (aZ) 4 , is not expected to improve the trend significantly, because i 

the ratio (d o/d ft) , . , _ to the Rutherford cross section still continues to 
third Born 

be higher than the exact phase-shift ratio to the Rutherford at small angles [8]. 
Figure 5 shows the results of the same calculations for an energy transfer of 
150 keV. 

Table 2 shows a comparison between the charge correction | Az | to 
the first Born approximation obtained by using the second Born approximation 
and the phase-shift calculations in equations (4) through (6). The energy 
transfer range used in equation (4) was assumed to be 50 to 150 keV. One notes 
from the table that corrections predicted by using the phase-shift formula range 
from about 1 for iron to 4 units of charge for Z = 75 and an incident nucleus 
with a kinetic energy of 1.457 GeV/nucleon. The second Born approximation, 
although agreeing up to Z » 52, predicts a correction of 10 units of charge at 
Z = 104. 


CONCLUSIONS 

Because the series used to calculate the phase-shift formula is slowly 
and conditionally convergent, one possible source of error would be that con- 
vergence has not been achieved, especially for those energy transfers that 
correspond to very small angles. However, great care has been exercised in 
the calculations to obtain convergence, both with respect to including transfor- 
mations to obtain a more rapidly converging series and by including a large num- 
ber (100) of terms. In addition, convergence tests were made in which up to 
250 terms were used . 

One notes immediately from the preceding results that the use of the 
second Born approximation, in the energy transfer region under consideration, 
leads to a gross overestimation in the correction to the Z 2 dependence of very 
highly charged nuclei producing knock-on electrons. 

One also finds the largest correction in charge determination due to the 
departure of the scattering cross section from a Z 2 dependence to be approxi- 
mately 4 units of charge at a kinetic energy of 1 .457 GeV/nucleon (/? = 0.92) . 

This is slightly higher than that quoted by Fowler [1] and much lower than the 
correction given by Semikoz [3]. 
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TABLE 1. CONVERGENCE TESTS FOR F^ G lt AND dcr/dT 
(T 0 is the incident nuclei kinetic energy in units of GeV/nucleon. 
da/dT is in units of Barns /Me V . ) 



a. Z = 92 - 0 = 0.90 (T 0 = 1.216) 

T = 50 keV 

T = 150 keV 

No Terms 

30 

100 

250 

30 

— 

100 

250 

Re F 4 

0.02434 

0.02381 

0.02381 

0.15859 

0.15509 

0.15509 

Im F t 

-0.16202 

-0.05584 

-0.05583 

0.02239 

0.01397 

0.01397 * 

Re Gj 

-5.40927 

-6.98770 

-6.98808 

-2.53794 

-2.53166 

-2.53166 

Im Gj 

sum 


-1.13001 

-3.88726 

-3.92990 

-3.92990 

dcr/dT 


1.1252x10® 



1.1252x10® 

1.3248x10 s 

1.3228x10® 

1.3228x10® 



b. Z = 92 - 0 = 0.97 (T 0 = 3.0) 

T = 50 keV 

T = 150 keV 

No Terms 

30 

100 

250 

30 

100 

250 

Re F t 

-2.62557 

-2.03047 

-2.02502 

0.36344 

-0.01404 

-0.01406 

Im Ft 

-28.50285 

0.60027 

0.59835 

-0.44164 

-0.03884 

-0.03883 

Re-Gj 

illlgg 


-12.82090 

-5.28019 

-7.96891 

-7.96992 

| 

■1 


BB 

-1.55627 

-1.29747 

-1.29659 





1.0380x10® 

1.0801x10® 

1.0802x10® 
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TABLE 2. COMPARISON OF |AZ| CHARGE CORRECTIONS TO THE FIRST 
BORN APPROXIMATION BASED UPON THE SECOND BORN APPROXIMATION 
(B2) AND MOTT PHASE-SHIFT CALCULATION (PS). 

(T 0 is the approximation incident nucleus kinetic energy in 
units of GeV/nucleon. Charge corrections shown are 
rounded to the nearest integer. ) 


\ z 
0 \ 

26 

52 

60 

75 

92 

104 

B2 

PS 

B2 

PS 

B2 

PS 

B2 

PS 

B2 

PS 

B2 

PS 

0.90 

(T 0 =1.212) 

1 

1 

3 

3 

4 

4 

B 

4 

B 

fl 

fl 

1 

0.92 

(T 0 =1.457) 

1 

1 

3 

3 

4 

3 

B 

4 

B 

B 

fl 

1 

0.95 

(T 0 =2. 068) 

1 

D 

3 

2 

3 

2 

5 

3 

B 

3 

8 

3 

IW 

- 

- 

2 

2 

2 

2 

4 

2 

5 

2 

■ 

2 
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da /dT (BARNS/MeV) 







Figure 2. Comparison between ratios of phase-shift calculation to 
(1) second Born Approximation and (2) first Born approx- 
imation^ /Incident nucleus has charge Z = 92 with /? = 0.92. ) 



Figure 3. Comparison between ratios of phase-shift calculation to 
(1) second Born approximation and (2) first Born approx- 
imation. (incident nucleus has charge Z = 26 with /3 = 0.92.) 
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Figure 4. Ratio of phase-shift cross section and Rutherford cross 
section as function of Z for 0 = 0.92 and /3 = 0.9712, with ratio 
of second Born approximation to Rutherford cross section 
included for reference. (The energy transfer is 50 keV. ) 



Figure 5. Ratio of phase-shift cross section and Rutherford cross 
section as function of Z for /3 = 0.92 and /? = 0.9712, with ratio 
of second Born approximation to Rutherford cross section 
included for reference. (The energy transfer is 150 keV. ) 
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